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Flicker (l/f) voltage noise spectrum is derived from finite-temperature quantum electromagnetic 
fluctuations produced by elementary charge carriers in external electric field. It is suggested that 
deviations of the frequency exponent 7 from unity, observed in thin metal films, can be attributed 
to quantum backreaction of the conducting medium on the fluctuating field of the charge carrier. 
This backreaction is described phenomenologically in terms of the effective momentum space dimen- 
sionality, D . Using the dimensional continuation technique, it is shown that the combined action of 
the photon heat bath and external field results in a 1//^ -contribution to the spectral density of the 
two-point correlation function of electromagnetic field. The frequency exponent is found to be equal 
to 1 + 5, where 5 = 3 — D is a reduction of the momentum space dimensionality. This result is applied 
to the case of a biased conducting sample, and a general expression for the voltage power spectrum 
is obtained which possesses all characteristic properties of observed fiicker noise spectra. The range 
of validity of this expression covers well the whole measured frequency band. Gauge independence 
of the power spectrum is proved. It is shown that the obtained results naturally resolve the problem 
of divergence of the total noise power. A detailed comparison with the experimental data on fiicker 
noise measurements in metal films is given. 
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I. INTRODUCTION 

As is well-known, power spectra of voltage fluctuations in all conducting media exhibit a universal low-frequency 
behavior ~ for sufflciently small frequencies / they scale as 1/ , with the frequency exponent 7 about unity. This 
asymptotic behavior is a manifestation of the so-called flicker noise present in any system containing free-like charged 
particle states [J. Although the value of 7 and an overall proportionality factor in this power law depend on 
many factors such as sample material and sample geometry, system temperature, etc., there are some important 
characteristic properties possessed by all flicker noise power spectra. Namely, it is well established experimentally that 
the noise produced by a biased sample is proportional to the voltage bias squared, and roughly inversely proportional 
to its volume. The enigmatic property of the l//-spectrum is its unboundedness. Flicker noise is present in the whole 
measured frequency band covering more than twelve decades. Experiments show no flattening of the spectrum for 
frequencies as low as 10~^ Hz. Although on the opposite side of the spectrum the l//-component is dominated by 
other types of noise (it compares with thermal noise usually at / ~ 1 Hz), it has been detected at frequencies as high 
as lO^Hz. 

Despite numerous models suggested since its discovery eighty years ago the origin of flicker noise still remains 
an open issue. There is a widespread opinion that this noise arises from resistance fluctuations, which is quite natural 
taking into account its dependence on the applied bias. It has been proposed that the resistance fluctuations possessing 
the other properties of flicker noise might result from temperature fluctuations [3|,yi , fluctuations of the charge carrier 
mobility 1, 10, Si, 

or of the number of charge carriers [l^. El!, El, El, Ej, llB, El|- AH these models, however, 
have restricted validity, because they involve one or another assumption specific to the problem under consideration. 
For instance, assuming that the resistance fluctuations arise from temperature fluctuations, one has to choose an 
appropriate spatial correlation of these fluctuations in order to obtain the desired proflle of the power spectrum. In 
addition to that, some models involve artificial normalization of the power spectrum, needed to come up with the 
observed noise level. 

Perhaps the main difficulty for theoretical explanation is the unboundedness of flicker noise spectrum. There are 
many physical mechanisms that generate noise whose power spectrum has the l//-proflle in some frequency domain. 
But these domains are so narrow in comparison with the whole measured band that the corresponding mechanisms 
cannot be considered as the general mechanisms of flicker noise generation. For instance, according to Ref. defect 
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motion in carbon conductors generates noise with power spectrum close to the inverse frequency dependence in the 
frequency range 10^ Hz to 10'' Hz, while outside this interval it switches to 1//^. At the same time, the omnipresence 
of flicker noise and high universality of its properties suggest that there must exist an equally universal reason for its 
occurrence. 

This source is naturally expected to have a quantum origin. Although some of the models suggested so far do 
consider various quantum effects as underlying mechanisms of flicker noise (such as, for instance, trapping of charge 
carriers), it may well be that its origin is to be sought at the most fundamental level. Namely, it is plausible that 
the phenomenon of flicker noise has its roots in the very quantum nature of interaction of elementary charges with 
electromagnetic field. From this point of view, the problem has been attacked by Handel [l^, who suggested that 
flicker noise is the result of low-energy photon emission accompanying any scattering process, and is related to the 
infrared divergence of the cross-section considered as a function of the energy loss. Later, the argument was modified 
and the so-called coherent quantum 1// effect described [l^, which is connected with the infrared properties of the 
dressed electron propagator. Although Handel's th eory has been severely criticized in many respects [20], it has found 
support in independent investigations of Refs. [2l|, |2^ . 

An essentially different quantum approach to the problem was proposed recently in Refs. ^23', '24'|. In this approach, 
flicker noise is treated as originating from quantum fluctuations of individual electric fields of charge carriers. As was 
shown in detail in Ref . [2^ , spectral density of the two-point correlation function of the Coulomb field produced by 
a charge carrier exhibits the above-mentioned characteristic properties of flicker noise. Namely, the low-frequency 
asymptotic of the spectral density is 1//, the noise intensity induced by external electric field is proportional to the 
field strength squared, and inversely proportional to the spatial separation between the field producing particle and 
the observation point. In application to the case of a conducting sample, the low-frequency asymptotic of the voltage 
fluctuation power spectrum is found to be 

Cuif)-^-, V=^9^^T, (1) 

where Uq is the voltage bias, T the system temperature, a the fine structure constant, fi the charge carrier mobility, 
c the speed of light in vacuum, and g a geometrical factor which is roughly inversely proportional to the sample size. 
The range of validity of this result turns out to be notably wide: The term "low-frequency asymptotic" means that for 
a given T, Eq. ([1]) is valid for frequencies satisfying / <C lO^^THz, with T expressed in °K. Covering well the whole 
band where flicker noise has ever been observed, this condition in particular sets no low-frequency cutoff, implying 
that the l//-spectrum extends down to zero. The fact that the found asymptotic does not require a low-frequency 
cutoff is the consequence of its oddness with respect to frequency. As discussed in Ref. [l^], appearance of such 
contributions to the spectral density is related to the inhomogeneity in time of fluctuations produced by individual 
charge carriers, and provides a natural resolution to the problem of divergence of the total noise power. 

It was demonstrated in Ref. [24] that Eq. ([1]) is in agreement with the experimental results of 1 //-noise measurements 
in metals. Naturally, in this comparison only genuine 1// noise data were used, i.e., the data that fit the law 1/ in 
which 7=1, within experimental error. Experiments show, however, that generally power spectra follow the one-over- 
f law only in sufhciently thick samples, while in thin samples (films, whiskers) large deviations of 7 (up to 7 = 1.5) 
are often observed. The purpose of this paper is to show that these deviations can be described within the developed 
theory by taking into account backrcaction of the conducting medium on the fluctuating electric field produced by a 
charged carrier. Staying within the one-particle picture of flicker noise generation, developed in Refs. (23l. |2^. this 
backrcaction can be described as an effective reduction of the momentum space dimensionality. It turns out that 
this reduction can be naturally realized using the well-known techniques of dimensional continuation [s^] . It will be 
shown that the power spectrum of electromagnetic fluctuations, continued in this way to 7 > 1, adequately describes 
the observed properties of flicker noise spectra in metal films. 

The paper is organized as follows. The effect of the backrcaction on photon propagation in a conducting film is 
discussed in Sec. Ill Al In Sec. IIIBl the influence of the heat bath on quantum propagation of electromagnetic and 
charged field quanta is considered, and the contributions relevant in the low-frequency regime are identified. The 
power spectral functions of the Coulomb field and voltage fluctuations are defined in Sec. Ill CI and written down 
in the form convenient for explicit calculations. Power spectrum of electromagnetic fluctuations in the presence of 
external electric field is evaluated in Sec. IIIII The low-frequency asymptotic of the voltage power spectrum is found 
to obey the l//'''-law with 7 = 3 — D, where D is the effective dimensionality of momentum space. Application of 
the obtained result to solids and comparison with experimental data is given in Sec. IIVI Some important auxiliary 
material used in the main text is collected in two appendices. Appendix A contains the proof of gauge independence 
of the voltage correlation function. Derivation of the dimensionally reduced photon propagator in the axial gauge is 
given in Appendix B. 

Units in which h — c — 1 are used throughout this paper except the end of Sec. IIIII and Sec. IIVI The spacetime 
metric 77^^ is defined mostly negative, i.e., 77^1/ = diag{-|-l, —1, —1, —1}. 
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II. PRELIMINARIES 
A. Photon propagation in a conducting film 

Consider a conducting film, i.e., a sample which is thin in one, say, 3-direction. The film will be assumed homo- 
geneous in (xi, X2)-plane, but otherwise arbitrary. Let the film be in a constant homogeneous electric field parallel 
to the (a:i, a;2)-plane, and denote U the voltage measured between two leads attached to the film at the distance L 
which is much larger than the film thickness a. We are interested in quantum properties of the electromagnetic field 
produced by a charge carrier moving in the film at finite temperature T. Specifically, finite-temperature correlations 
in the values of the particle's Coulomb field in the presence of external electric field will be investigated. As long as 
backreaction of the conducting medium on the fluctuating field of the charged particle is neglected, these correlations 
are a single-particle effect, in the sense that to the leading order in the electromagnetic coupling, only fields produced 
by one and the same particle correlate. However, this is no longer the case upon account of the backreaction. The 
point is that under usual conditions of flicker noise measurements, the backreaction is to be considered as an effect of 
zero order in the electromagnetic coupling. Indeed, since the characteristic time of the charge density response to the 
field fluctuation is much smaller (normally, by ten orders of magnitude at least) than the time 1// of the field mea- 
surement, this response can be considered quasi-stationary, leading to one and the same charge density redistribution 
independently of the value of charge carried by elementary medium constituents. Since these redistributions tend to 
compensate the field fiuctuation, the backreaction leads to damping of the electric field fluctuation. Concerning their 
effect on quantum propagation of the particle's field, significance of these redistributions depends on the mechanism 
governing the response. From this standpoint, response of the medium through the ordinary electric conduction is 
inconsequential because this is a classical effect, in the sense that it relates quantities averaged over many decoherent 
particles - the mean electric field and mean electric charge density. Things differ, however, upon account of specifically 
quantum mechanisms such as the exchange interaction between charge carriers. Since this interaction is independent 
of the particle charge, it changes the structure of electromagnetic correlations already at the lowest order. This inter- 
action leads to correlations of changes in the charge density distribution, induced by the fluctuating electric fleld of 
the given particle. Being related to the symmetry properties of the charged particles state, the exchange interaction 
leads to density correlations which are coherent, and therefore so is the response of the medium on the fluctuating 
field. It is important that the charge carriers correlated this way react coherently on each single photon propagating 
in the medium, thus changing quantum properties of this propagation. In other words, the exchange charge density 
correlations modify the form of the electromagnetic field propagator. 

For definiteness, in the rest of this section we consider metal conductors. Since electronic component in metals 
is degenerate at ordinary temperatures, for a rough estimate of the exchange effect it will be sufficient to use the 
simplest model of ideal degenerate Fermi-gas of neutral particles. It is known [2^ that the exchange correlations of 
density fiuctuations in this case are destroyed by the thermal effects at distances 

_ hpF 
2TrmT ' 

where is the Fermi momentum, and m the particle mass (for distances r > tq, density correlations fall off 
exponentially as e~'"/''°). For electrons in a metal, pp can be estimated as h/d, d ^ 10~^ cm being the lattice spacing, 
and hence 

To « 10~^ cm . 

At shorter distances, the correlation function of density fluctuations scales as r~'^, and at r « d fluctuations become 
completely correlated. Thus, one can expect that the propagation of photons polarized in the 3-direction will be 
partially suppressed by the exchange effects in films with thickness a ro, and damped almost completely in the 
case a « d. This applies to the real transversely polarized photons as well as to the longitudinal "photons" describing 
Coulomb fields of charged particles. Let us turn to the consequences of this damping regarding the form of the 
electromagnetic field propagator, 

Df,,{x ~ x') ^ i{7Af,ix)A,{x')) , 

where T denotes the usual time ordering of field operators, and averaging is over the given photon state (equilibrium 
distribution at temperature T). The quantum suppression of the 3-component of the electric field means that the 
electromagnetic field operator is constrained by the condition 

E 1 OAs dAo _ ^ 
^ c dt dx3 
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At the same time, it is well known that the form of the photon propagator is to a certain extent arbitrary. This 
arbitrariness reflects the freedom in choosing gauge conditions used to fix the gradient invariance. On the other hand, 
it is proved in Appendix A that the voltage correlation function is gauge- independent, so the choice of the gauge is 
at our disposal. In the present case, it is convenient to choose the axial gauge 



Then the above operator constraint simplifies to 



^3-0. (2) 



§^2^0, (3) 

0X3 



This condition means that the normal mode decomposition of the scalar potential does not contain wave vectors with 
nonzero 3-component. Thus, in this particular gauge, the suppression of E3 implies that the "temporal" photons do 
not propagate in the 3-direction. In other words, the momentum space propagator of these photons becomes effectively 
two-dimensional. It is important, furthermore, that the complete photon propagator turns out to be diagonal in the 
axial gauge under the condition ([3]). Namely, the zero-temperature photon propagator is found in Appendix B to 
have the following form 

- iv,. + n,n^) \ I , (4) 

where is the unit vector in the 3-direction. As a consequence, the finite-temperature propagator is also diagonal 
(see the next section). In view of this fact, Dqq turns out to be the only component of the photon propagator that 
determines the power spectrum of voltage fluctuations produced by nonrelativistic charge carriers. This is because 
whenever D^i, is contracted with the electromagnetic current J'^ , contribution of the components with v ^ is 
suppressed by the factor \q\/mc, where q is the charged particle momentum, so that D^^J'^ sa D^qJ'^ = {DooJ^, 0, 0, 0). 
In what follows, therefore, we will be dealing only with this component of the photon propagator, and drop the lower 
Lorentz indices, for brevity. 

Our choice of units h — c — 1 still leaves freedom in choosing the unit of length. It proves to be convenient to fix 
it by setting a = 1, which will be assumed until the end of Sec. IIIII 



B. Finite-temperature contribution to the propagators 



Let us now consider the influence of the heat bath on the photon and charged particle propagation. As long as 
the mean electric potential produced by an elementary charged particle is considered, the photon heat bath has no 
effect. This is because the 4-vector of momentum transfer to a free massive particle, p, is always spacelike, < 0, so 
the distribution of real photons appearing in the definition of the photon propagator is irrelevant. The consequence 
of this is that the quantities built from the mean field, such as the disconnected part of the correlation function, are 
not affected by the photon heat bath. Things change, however, when the connected part of the two-point correlation 
function of electric potential is considered. It is defined by the following symmetric expression 

CoT(^;^') - l{in\Ao{x)Ao{x') + Ao{x')Ao{x)\in) , (5) 

where x and x' are the spacetime coordinates of two observation points, Aq is the scalar potential operator in the 
Heisenberg picture, and |in) denotes the given in state of the system "charged particle -I- electromagnetic field." In 
the two-photon processes, photon momenta are allowed to take on lightlike directions, and hence the photon heat 
bath does contribute to the function Cqq'^{x] x'). 

As is well-known, the ordinary Feynman rules of the S- matrix theory are not generally applicable for the calculation 
of in-in expectation values, and must be modified, e.g. , according to Schwinger and Keldysh |26l . \2T\ . This complication 
was overcome in (23 | by rewriting Eq. ([5]) in the form 

CoT(^; = Re(out|r{io(a;)io(:r')}|in) , (6) 

which allows the use of the S-matrix rules. This transformation uses equivalence of the one-particle in and out states 
(and the Hermiticity of the electromagnetic field operator). Thus, in order to calculate the in-in expectation value 
^ taking into account the heat bath effect, the standard finite-temperature-field-theory techniques can be used [28l |. 
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Below we employ a version of the real time formalism, developed in [29|, which is especially convenient in actual 
calculations since the momentum space propagators in this formulation do not involve the step function. 

The real time formulation involves doubling of all fields, which will be specified by a two-valued lower index. 
According to the diagrammatic rules derived in 29], the photon propagator at finite temperature can be obtained 
in momentum space (of arbitrary dimensionality) by replacing the quantity 47r/(fc^ + iO) in Eq. ([4]) by the following 
matrix 

^(k)-An( ^''^^^ ^''^^^] (7) 



where 

/?ii(fc) = -D*^^{k) 



1 27ri(5(fc2) 



= i^2i(fc) = - '"jj_l ' 

In applications to the problem of l//-noise considered below, the value of the product /3|fco| turns out to be very 
small. For instance, even for frequencies as large as lO^Hz, and temperatures as small as 1°K, it does not exceed 
h ■ 10^/fc K, 10~^ (fc is the Boltzmann constant), so the denominators in the above expressions can be replaced by 
/3|fco|, implying that the second term in Dn dominates. On the other hand, the temperature effect on the propagation 
of massive particles is much less prominent. For instance, in the case of conduction electrons in a crystal (this case 
will be used throughout as a standard example), the particle energy is of the order {h/dY /m. Setting m « 10~2^g, 
we find that {h/dY /mT ~ 10+^/T, with T expressed in °K. Hence, the temperature contribution can be completely 
neglected (for fermions as well as for bosons), and the propagator taken in the simple diagonal form 





Dt2{k) 



(8) 



DUk) = -Dt;{k) = {m^-k^~iQ)-\ 
This is for a scalar particle described by the action 

= [d^x{{df,^* +ieAf,(j)*)id''(l)~ieA''^)-m^(l)*(l)} , (9) 



where e is the particle charge. Account of particle spin, though adds some extra algebra, does not change the long- 
range properties of its field, so that following Ref. j24| we work with the simplest case of zero-spin particles. The 
matrix propagators are multiplied in the interaction vertices, generated by the triple and higher order terms in the 
Lagrangian, with an additional minus sign for the product of 2-components, as in the Schwinger-Keldysh techniques. 
The arguments of the Green functions are treated as 1-component fields. Finally, external particle lines represent 
normalized particle amplitudes or their conjugates, according to whether the particle is incoming or outgoing, just 
like in the conventional techniques. 

As was mentioned above, the photon propagator is dominated by the temperature contribution (as long as the 
range of momentum integration contains lightlike directions, see discussion in Sec. while in the massive particle 
propagator this contribution is negligible. It is important, on the other hand, that the heat bath affects significantly the 
real particle propagation, i.e., external matter lines in the diagrams. Bilinears of the particle amplitudes representing 
these lines are expressed eventually via statistical distribution function (see Sec. Illll for details). Thus, apart from 
explicit T-dependence coming from the photon propagator, the correlation function also depends on temperature 
implicitly through the particle statistical distribution. 

C. Power spectral densities of potential and voltage fluctuations 

Connected contribution to the power spectral density of electric potential fluctuations is obtained by Fourier trans- 
forming Eq. ^ with respect to the difference of the time instants t, t': 

C{x,x',t',Lu) = [ dTC5°"(a;,i' + r;a;',i')e"'"^. (10) 
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The upper and lower indices in the notation of the correlation function are suppressed in the left hand side, for brevity. 
We are interested ultimately in the power spectrum of voltage fluctuations, C[/, measured between two observation 
points x,x' (the two leads attached to the film). The connected contribution to the voltage correlation function is 
given by 

Cu{x,x',t,t') = ]^l:in\U{t)U{t') + U{t')U{t)\m) , (11) 

where U{t) — Ao(a;,f) — Aq{x' ,t) is the operator of voltage between the two points. This function is separately 
symmetric with respect to the interchanges x <-> x', and t <-!■ t' , unlike the function Cqq"(x;x') which is only 
symmetric under x *^ x' . Substituting the definition of U{t) in Eq. (Ilip . the former can be expressed via the latter as 

Cjj{x,x,t,t) — Cqq'^ {x,t; x,t ) + Cqq'^ {x ,t] X ,t ) 

- [C-oo''i^,t- x', t') + Qo'^i^', t;x,t')] . (12) 

Accordingly, the power spectrum of voltage fluctuations, defined by 

Cu{x,x',t',u;)= J dTCu{x,t' + T,x',t')e-'"-' , (13) 

— oo 

is expressed through that of potential fiuctuations as 

Cu{x, x', t', uj) = C{x, X, t', uj) + C{x\ x\ t', uj) - \C[x, x' , t', lu) + C{x' , x, t' , uj)] . (14) 

Although Cu{x,x' ,t,t') is symmetric with respect to the interchange t <-> t' , it depends on both time arguments 
separately, and therefore Cij{x^ x' , t,uj) does not have to be an even function of w. 

When calculating the power spectrum of potential fluctuations according to Eqs. ©, (fTU|) . it is convenient to 
perform the Fourier transformation under the sign "Re" in Eq. ([6]). For this purpose, we introduce Fourier transform 
of the two-point Green function: 

G{x,x',t',uj) = [ dTGix,t' + T;x',t')e-''^^ , Gix; x') = {out\T{Aoix)Aoix')}\m) , 



(15) 

with the help of which the power spectrum of potential fluctuations can be written as 

C{x,x',t',uj) = ^Re{G{x,x',t',uj) + G{x,x',t',-uj)} 

+ ^lm{G{x,x',t',Lu)-G{x,x',t',-uj)} . (16) 

We see that contributions to the function C{x,x' ,t' ,uj), and hence to the voltage power spectrum, can be either real 
even, or imaginary odd functions of frequency. 



III. LOW-FREQUENCY ASYMPTOTIC OF THE POWER SPECTRUM 

In this section, spectral density of the voltage correlation function will be evaluated in the low-frequency limit 
taking into account the flnitc-temperature effect and the influence of constant homogeneous external electric field. 
Let us first discuss the role of the particle collisions in this calculation. It was mentioned in Sec. Ill Al that, barring 
the exchange interaction, correlations in the values of the electromagnetic fields produced by the charge carriers is a 
one-particle effect within the leading order in the electromagnetic coupling, in that only fields produced by one and the 
same particle correlate. At the same time, in the presence of the external electric field the direct particle interactions 
cannot be neglected completely. Indeed, although this field leads to a relatively small corrections to the charge carrier 
wave function and its propagator in all practically relevant cases, the effect of constant homogeneous field on the 
free-like particles cannot be treated perturbatively. The role of the particle collisions is to prevent the charge carrier 
from gaining too much momentum from the field, thus cutting down its effect. This is pictured schematically in Fig.[T] 
where the particle collisions are symbolized by a virtual photon interchange between the particles. The external field 
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and the collisions affect both the particle wave function and its propagator. As to the former, account of these two 
factors is accomplished by replacing the particle momentum probability distribution by the statistical distribution 
function, obtained as a solution of the kinetic equation in the presence of external electric field (this point will be 
discussed in more detail later in this section). Regarding the particle propagator, however, the issue is not that simple, 
because in order to extract the low-frequency asymptotic of the power spectrum one needs an explicit expression for 
the propagator, incorporating both effects, which is unknown. This difficulty can be overcome by the following trick. 
As was mentioned above, the external field gives rise only to a relatively small correction to the particle propagator. 
It is, therefore, sufficient to determine this correction in the linear approximation, i.e., to the first order with respect 
to the electromagnetic coupling. We now make an assumption (confirmed by the result of the calculation) that the 
low-frequency asymptotic of the power spectrum is independent of the charge carrier mass. Then this asymptotic can 
be found formally in the large mass limit. For a given frequency uj, one can always take the mass large enough so 
as to justify perturbative treatment of the external field effect. To the first order in the electromagnetic coupling, 
corrections to the particle propagator due to external electric field and particle collisions are well-separated from each 
other because they are superimposed linearly. On the other hand, since the characteristic time of particle collisions is 
very small, their contribution to the particle propagator is inconsequential in calculating the low-frequency asymptotic 
of the power spectrum. Thus, the lowest order contribution to the correlation function is represented by diagrams 
with a single insertion of the vertex describing interaction of the charged particle with external electric field. These 
diagrams are shown in Fig. O Distinguishing the contributions of the diagrams [21 a) and [D(b) by the corresponding 
Latin subscript, we have 

Ga{x,x') = -e^ /// d^zd^z'd^z"Ao{z") 

\{z'-x')\ , (17) 
J 11 



X i^{x - z) 
Gh{x,x') = Ga{,x' ,x) 

where 



^{z) do S^(z - z") < 2)^(z" - z') ^P*{z') 



^ I d 

ipda-tjj = ipdoTp - i^doifi , = — , 

CZq 

Ao(z") = —{Ez") is the external field potential,^ and the given particle state. Next, we go over to momentum 
space with the help of Eqs. ([7]), ([8]), and introduce the spectral function for G{x,x') according to Eq. (|15p . Assuming 
also that the charged particle is nonrelativistic, and taking into account that — 0, we write the matrix product 
longhand^ 

xe'P°(*'-*«)-Hp-'=')^'ylo(fc')>/a(p, q, k\ x - x\ uj) , (18) 

G,{x,x',t',u:)^G:{x,x',t',-cu), ^ fc^l! - |1 , (19) 

2m Zm 



where 



Ja{p, q, k', x-x',uj)^- ie{2mf J A^e'M-'-) (20) 

X ^^Dii{k)Df^{q + k)Df^{q + k + k')Dii{k + k'^p) 
^Di2{k)Dt2{q + k)D^2iq + k + k')D2i{k + fc' -p)} 
Ao{k')^-t{27r)'^(^E^y('^\k'), k' ^(Q,k[,k'2,0), D = 2. 



kn—u 



^ We do not include an arbitrary constant in this expression because it is clear in advance that it cannot affect the final result. It is 
proved in Appendix A that the correlation function is actually invariant under the most general gauge variations of the electromagnetic 
potential. 

^ The relation Gf,{x, x' ,t' ,Ld) = G*{x, x' ,t' , —lo) is conveniently proved using the following sequence of substitutions in Eqs. HISII - II20I I: 
fe— >fc-|-p — fc',q— >qr — p, and then p —p. The extra factor (—1) coming from the complex conjugation of the imaginary unit in 
Ao{k') is compensated by that from the integration by parts with respect to fc'. 
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Here is the charged particle 4- momentum, and ij{q) its momentum wave function at some time instant t^^ normal- 
ized by 

"^'^1^(9)1^ = 1. (21) 



That the integrations over momenta k,p,k' in the expressions (US]), (PO)) turned out to be two-dimensional follows 
from the expression for the photon propagator, derived in the preceding section, which in the current notation can be 
written as 

= j (^^rp+T^^^'^sw, © = 2, 

S) being given by Eq. ([7]). As we saw in Sec. Ill Al this reduction of momentum space dimensionality from 2) = 3 to 
2) = 2 is the consequence of quantum damping of temporal photon propagation in the 3-direction. This suppression 
is noticeable only in sufRciently thin films with a < ro, being formally complete in the practically unattainable case 
a = d, and is negligible in thick samples with a 3> tq. To interpolate between these two limiting cases, it is natural 
to consider D in the above formulas as a phenomenological parameter allowed to take on non-integral values between 
2 and 3, and which therefore can be called "effective dimensionality" of the film. The difference 3 — D = (5 thus will 
be an effective reduction of the film dimensionality, caused by the suppression. To accomplish this interpolation, we 
employ the well-known technique of dimensional continuation, developed by Wilson (in condensed matter theory), 
and by t'Hooft and Veltman (in quantum field theory) 'so']. 

Before evaluating the D-integrals, let us note that the integrand in Eq. ([^0]) can be considerably simplified. First of 
all, the second term in the curly brackets can be neglected. Indeed, in view of the factor Di2(k) which is proportional 
to (5(fc^), and the condition fco = the momentum k contributes only a tiny value to the argument of the factor D21] 
for electrons in a crystal, for instance, the ratio |fc|/|p| ~ {hw / c) / {h/ L) = ujL/c w lO^^^wL, where L are expressed 
in the CGS system of units (L is the sample length) . Even for L as large as 1 cm, this ratio is very small for all 
practically relevant frequencies. Taking into account also that the momentum k' is set eventually equal to zero, the 
factor D2i{k + k' — p) can be written simply as I?2i(p) ^ ^{p^)- But the argument of this delta-function is always 
nonzero, because momentum transfer p to the massive particle is spacelike. Furthermore, using explicit expression for 
the photon propagators, their product in the first term in the curly brackets reads 

1 1 27ri(5(fc2) 1 

7^ii(fc)Dn(fc + fc'-p) = ^ ' 



[fc2 + iO] [{k + k' - p)2 + iO] e'3|'=ol - 1 (fc + /c' - p)2 + io 

1 2TTiS{{k + k' ^ pf) 2TriS{k^) 2m5{{k ^ k' - pf) 

k^ + iO g0\ko-po\ _ I gPlkol _ I g/3|fco-pol _ I 

As before, the last term in this expression can be omitted, while the first term is dominated by the second, as we 
saw in Sec. IIIBI Furthermore, kg enters the temperature exponent in the third term in the combination (fco — po), 
and therefore, this term does not contribute to the leading term of the low-frequency asymptotic, because in practice 
\'^\ ^ Po- Indeed, estimating the energy transfer po {pq)/m ~ h^/mdL, one finds for our standard example 
I'^/t'oI ^ 10~^\uj\L. Finally, in the product of the second term with Df^{q + k)Df^{q + k + k'), the pole of the function 
-^11(9 + ^) does not contribute because of the factor S{k'^). This is again a consequence of the requirement that 
the momentum transfer to the massive particle on-shell be spacelike: conditions k^ = 0, q^ = m? , and {q + fc)^ = 
m? cannot be satisfied altogether. Hence, the scalar particle propagator in this product can be written simply as 
Dti{q + k) = -l/2{qk) « -l/2muj. Setting also (e'^l'^l - 1)-^ « T/\uj\, one thus finds 



Ga{x,x',t',uJ)^-^{2n)^{4ne)' /// ^^^^(g)V/(g + p) 



xe ■ 



d^q d'^p d'^k' 
(2^)3 (27r)^ (2^)^ 
d_ 
dk' 



where 



47remr f d^k pf^j^g + k + k')S{k^) ^ 

ko—LU 



Mp,q,k,x~x,u;) = ~—^J^e ^k + k'~pr + ^0 



The singular at w = contribution to the function Ga{x,x' ,t' ,uj) comes from integration over small k. The strength 
of this singularity is determined by the poles of the propagators in the integrand, and we have to decide which of 
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them gives rise to the strongest singularity after performing the fe'-integration. Consider first the case when the 
fc'-derivative acts on the factors [{k + k' — p)^ + iO]~^, e^^^P^'' •''^ . Since these depend on the difference (p — k'), 
changing d/dk' — > —d/dp, and then integrating by parts with respect to p in Eq. ((22|) . this derivative is rendered to 
act on terms independent of fe. Thus, of all the fe'-dependent factors in the integrand only Df^{q + k + k') is to be 
differentiated in effect. This brings the above expression for Ga to the form 



(27r)3 (27r) 
xe*P°(*'-*o)-*P-' {Eq) Ja{p, q,x- x' , cu) 



with 



ft is seen from Eq. that if a term in the function C{x, x' , t' , oj) is independent of one of the arguments x, x', then 
it does not contribute to the voltage power spectrum. Hence, we expand the exponent in the integrand of (P^ . and 
retain only one term leading in the low-frequency limit. Using the rules of ©-integration ^3d\, one finds 

Mp,q,x-x\u;) = ^£^1 ^[kix'~x)f6i^'-k') 

■weTix' -xf f (Pk fc2 

' -5{uj - k^) 



where 



TTeT{x'-xf |a;p5D 



r('D/2) 

is the area of unit hypersphere in D dimensional space (r(x) is the Eulcr function). Thus, 



■Dr(H/2) 2771^^1^^1-5 J J (27r)3 (27r)2'^'^'^ p2 ' 

where L — \x' — x\. The total contribution to the spectral density of the two-point Green function is 

G{x, x\ t', oj) = Ga{x, x', t', oj) + G* (cc, x', t', -uj) . 

Substituting this into Eq. (jl6p one sees that since Ga{x, a;', lo) is odd with respect to frequency, it is the imaginary 
part of the Green function that contributes to the correlation function C(cc, a;', t', uj) which therefore takes the form : 



DrCD/2) mLo\uo\^ J J (27r)3 (27r)25 ^^^^^ \h ' i^j^ ^2 ■ 

(22) 

In a many-particle system, this result is to be expressed through the one-particle density matrix, which is accomplished 
by replacing ip*{q')4'iQ) ~^ 0a{Q',Q), where go is the momentum space density matrix at the time instant Iq. Recall 
that to is the instant at which the particle state ipiQ) is prepared. It can be identified, for instance, as the moment 
the charge carrier enters the sample, or escapes from a surface trap, etc. The factor e'^"^* in the integrand 
realizes evolution of the density matrix from the instant to to t'. Since Po = {q + p)^/2m — q^/2m, the product 
gipo(t ~*o)go(q +p,q) describes a particle evolving freely on the interval {to,t'). On the other hand, as was already 
mentioned in the beginning of this section, in order to justify perturbative treatment of the external field effect on 
the real particle states, it is necessary to take into account particle collisions. For this purpose, it is sufficient to 
consider these collisions as instantaneous. Then the interval {to, t') is divided into a sequence of short time intervals of 
duration Tf (the particle mean free time) , on each of which the density matrix evolves freely, and changes abruptly at 
the collision instants. Going through this sequence, the density matrix tends to the stationary statistical distribution 
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function, q), which is independent of the initial particle state. What is important here is the sign of the difference 
(i' — to)- Recall that t' is a fixed time instant to count off the time interval r with respect to which the correlation 
function is Fourier-transformed, and that each particle has its own t^. This means that for a given the system is 
observed during the time interval (t' — At, t' + Ai), where At ^ 1/t^, and t^s are distributed uniformly over this interval. 
The density matrix evolves forward (backward) in time, if t' > tg {t' < to)- But time reversal involves inversion of 
particle momentum, and therefore, the reciprocal contributions to the function G{x,x' ,t' ,lu) have opposite signs. To 
be more specific, let (t' — to) > 0. Then the exponent e*^°*^* realizes forward evolution of the density matrix, so 
that the integral in Eq. (|^^ takes eventually the form 

On the other hand, if (t' — tg) < 0, then the density matrix evolves backward. In momentum space, the initial state 
of the reversed motion is represented by the amplitude ^'{q) — ^*{~q). Taking complex conjugate of the integral in 
Eq. ([22|) (which does not change the value of G in view of the sign "Re"), and changing the integration variables 
q —q, p — > —p gives in this case 

'''' '''^,V3(9)^*(9 + P)^e-(*0-0-<P.'. 



(27r)3 {27r)^' '^'^ p 

Replacing ip*{q + p)ip{q) g{q + p, q), where g plays the role of momentum density matrix at the moment t', the 
exponent e^Poito-t ) governs forward evolution of this state on the interval {t',to), so that the above expression takes 
the form 

d^q (Pp . AEq) 



giq + p, g)^^e-'f" 



(2^)3(2^)!^'^^^ ' ^'^^ p2 

The density matrix here is the same as in (P^ . because the statistical distribution is independent of the initial state. 
We see that reciprocal contributions to the function G{x,x' ,t' ,uj) cancel each other when summed over all particles 
in the system. Thus, we arrive at the important conclusion that the total noise intensity is independent of the number 
of particles, and remains at the level of individual contribution. As was shown in Ref. [24j, this conclusion is also true 
of the disconnected part of the correlation function, though by virtue of quite different reasons. 

It is customary to further express the function g{q',q) via the real mixed distribution function, n{r,q), according 

to 



Q{q + p,q) = J 



Probability distributions for the particle position in a sample or its momentum can be obtained by integrating n{r,p) 
over all p or the sample volume, respectively. Using this in the expression (j23p . and substituting the latter into 
Eq. ^ yields 



■Dr(H/2) muj\uj\s JJJ (27r)3(27r)^ p^ V ' 2 

Shifting here q q — p/2, and dropping the purely imaginary term proportional to (Ep), the triple integral becomes 
purely real, so the symbol "Re" can be omitted. Integrating then over p with the help of the formula 

/ rf^P,.(..)4._ r(D/2-l) 



we obtain 



, , , 2i(27r)'^ e^TL"^ f f d^q n ,^ n(r,q) 
C{x,x\t\uj) ^ ^.^ --^ // —-^Sr{Eq) ^ 



/|D-2 



D{T> - 2) mu)\oj\^ J J (27r)3 ^ \r ~ x- 

Finally, substitution of this expression into Eq. p4p gives low-frequency asymptotic of the power spectrum of voltage 
fluctuations 

c„(...-,,-..) . - ^g:';, j^MEMr)) (^^ + , (24) 
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where 



is the local drift velocity of charge carriers, Q. denoting the sample volume. For a crystal in a homogeneous external 
field, U is a function of the crystalline direction, 

Vi = ^iikEk, i,fc = l,2,3, 

where is the charge carrier mobility tensor. To write down the final expression for the voltage power spectrum, 
we restore the ordinary units. Then Eq. (|24p takes the form, for lo — 27r/ > 0, 



where 



C^(a.,.',t',^) = -z-JL, .^-^^^,Tg^ (25) 



E 

= fJ-ikn-iUk , n = — , 
\^\ 

Ua ~ \E\L is the voltage bias applied to the sample (it is assumed that E \\ (a; — x'), as is usually the case in practice), 
and g a geometrical factor 

' ^ d'r ( j-^-^ + j-^—,) . (26) 



For w < 0, the sign of the right hand side in Eq. (|25p is opposite. If Fourier transformation is defined in a purely real 
form, i.e., as a decomposition in cos(u;r), sin(ci;r), rather than in e"^'^, then the spectral density is also real: 



Cu{x,x',t',u;) = j^. (27) 



It is seen from Eq. (|26p that the g-factor has a pole at D = 2, which means that the noise amplitude is unbounded 
in the limit 2) ^ 2. It was already mentioned in Sec. IIIII that the case D = 2 would correspond to the practically 
unreachable film thickness a — d. We now see that this case actually cannot be realized even theoretically. It should 
be emphasized in this connection that since in practice the film thickness always largely exceeds the lattice spacing, 
only a relatively small fraction of charge carriers in the film is correlated by the exchange interaction. Therefore, 
by continuity, the effective film dimensionality must be close to 3. In other words, the range of applicability of the 
developed theory is in any case limited to fiicker noise spectra characterized by sufficiently small values of (5 = 7 — 1. 

We mention for future reference that if the sample is an elongated (say, in x-direction) parallelepiped with the leads 
attached to its ends, then the g-factor can be evaluated approximately as 

\adx 2 ^28) 





where w is the sample width, and it is assumed that a < w ^ L. This is for d > 0. For d close to zero, this formula 
is not applicable, because the x-integral diverges near the sample ends. Cutting off the integral a,t x ~ w, one finds 
with logarithmic accuracy 

L 

wadx 2 L „ , ^ 

In - , ^ w . (29) 



3Lwa I X 3L w 



We note also that in the CGS system of units, the >f-factor reads 



1.62.101"-4*^^-22^mT, (30) 



where the absolute temperature T is to be expressed in °K. 
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IV. APPLICATIONS 
A. Unboundedness of flicker noise spectrum 

In this section, a special feature of the derived expression for the power spectrum, namely, its oddness in frequency, 
will be discussed in connection with the problem of observed absence of frequency limits of the l//'''-law. As was 
mentioned in Introduction, flicker noise has been detected in a very wide frequency band ^ 10^^ Hz to 10^ Hz. This 
fact represents one of the essential difhculties for theoretical explanation, because all physical mechanisms underlying 
existing models of flicker noise work in much narrower subbands, and none of the models suggested so far has been 
able to explain the observed plenum of the l//-spectrum. 

On the other hand, existence of bounds on this spectrum is generally believed to be necessary in order to guarantee 
finiteness of the total noise power. There is a well-known argument [3i| according to which these limits are actually 
unnecessary when the flicker noise exponent 7 is strictly equal to unity, because the logarithmic divergence of the 
total power is not a problem in this case in view of the existence of natural frequency cutoffs such as the inverse 
Planck time and lifetime of Universe. However, this reasoning does not work for 7 ^ 1, in which case divergence is a 
power of the cutoff. At the same time, the results obtained above reconcile unboundedness of l//-spectrum with the 
requirements of stationarity and finiteness of the total noise power in a quite natural way. Indeed, using Eq. (j25p we 
find that for S < 1 the integral 

+OC +OC +00 

duj „ , ,^ f doj „ , , , , , ^ n f sin(27r fr) 

— Cu{x,x',t',u:)e'^^ = 2i j —Cu{x,x',t',u;)sm{iUT)^2;<USj d/— 

-00 

converges in both limits / — > and / — s- 00. In particular, the singular contribution to the voltage variance (i.e. to 
the quantity Cu\t=o) vanishes. It is worth to note also that the total noise power would diverge for 6=1 even if >c 
were bounded at D = 2. 

Since appearance of odd contributions to the power spectrum is somewhat unusual in macroscopic fluctuation theory, 
let us discuss it in more detail. Under stationary external conditions, the voltage noise power spectrum (to be denoted 
below simply as Cu{t,t'), with the spatial arguments suppressed, for brevity) must be independent of t'. This is an 
expression of the noise stationarity, or, using a term more suitable for the subsequent discussion, time homogeneity 
with respect to the macroscopic system. It is usually realized as the requirement that Cu{t,t') be a function of the 
difference t — t' = T. Since Cu{t, t') is also symmetric with respect to the interchange t t', an immediate consequence 
of this is that it is actually a function of |r|, and hence the spectral density is a real even function of frequency. It is 
important, on the other hand, that time homogeneity is not necessarily exhibited by individual contributions to the 
total voltage fluctuation, whatever mechanism of flicker noise generation be. In particular, this property evidently 
does not take place at the microscopic level, i.e., with respect to elementary processes such as charge carrier trapping, 
surface or grain boundary scattering, etc. Stationarity of the macroscopic process emerges usually upon summation 
over a large number of individual contributions, so that this microscopic inhomogeneity turns out to be inconsequential. 
However, this summation is not the only way to obtain a stationary correlation function symmetric in t, t' . Another 
possibility, which is realized in the present paper, is that flicker noise may be a one-particle phenomenon, in the sense 
that the entire effect can be ascribed to elementary fluctuations produced by single charge carriers. In this case the 
function Cu{t, t') does not have to depend solely on |t|, and as the explicit calculations of Sec.lIIIIshow. it actually does 
not. As was mentioned above, elementary processes are inhomogeneous in time, and hence the symmetry with respect 
io t t' imposes no restriction on the r-dependence of the correlation function. The only remaining requirement, 
namely reality of the correlation function, implies that contributions to the spectral density must be real even, or 
imaginary odd functions of frequency [Cf. Eq. (jl6p ]. These two cases correspond to the Fourier decomposition of 
the function Cu(t' + r, t') in cos(ti'T) and sin(cL;T), respectively, and describe the parts symmetric and antisymmetric 
with respect to the difference of its time arguments. Finally, transition to the statistical distribution removes the 
t'-dependence of the power spectrum [Cf. discussion after Eq. (|22p ]. This restores macroscopic time homogeneity of 
the correlation function, but leaves the possibility of being odd with respect to the difference of its time arguments. 
In other words, dependence of the power spectrum on t' shows itself only at microscopic scales, while macroscopically 
fluctuations look as if they were homogeneous in time. 

The l//-spectrum derived in the previous section has no lower frequency cutoff. As to the upper bound, it is given 
by the condition / <C T [see Sec. IIIBj . or in the ordinary units, / <€. kT/h « lO^^THz, with T expressed in °K. We 
see that from the practical point of view, the obtained spectrum has no upper cutoff cither. 
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B. Comparison with experimental data 

Let us continue verification of the obtained result and show that Eq. (|25[) is in agreement with the other exper- 
imentaUy estabhshed properties of flicker noise. We will first discuss some general qualitative properties of flicker 
power spectra, predicted by Eq. (j25p . and then give a detailed quantitative comparison of these predictions with 
experimental data. 

1. Qualitative comparison with the experiment 

First of all, the power spectrum of quantum electromagnetic fluctuations, given by Eq. (|25p . is quadratic in the 
applied bias. This is perhaps the most solidly established property of flicker noise. Second, the noise level is generally 
inversely proportional to the sample size. Namely, the (7-factor describing dependence of the noise intensity on the 
sample dimensions increases roughly as a power of decreasing sample length or thickness, the exponent depending 
on the sample geometry as well as on the effective sample dimensionality D (which itself depends on the sample 
thickness.) As to the dependence of flicker noise amplitude on sample dimensions, agreement in the literature is not 
that good. Experiments are usually arranged so as to prove one of the two main competing points of view on the 
flicker noise origin, namely wether it is a bulk or surface effect. Although this issue is far from being resolved, there is 
no doubt that the noise level increases with decreasing sample size. To be more specific, we note that if 5 = const, then 
Eqs. (pS)) tell us that the noise produced by an elongated sample is proportional to a^^^L^^^. For experimental 
verification of this prediction we refer to [s^ , [33| which report the results of flicker noise measurements in various 
copper films with o = 400 A-;-2000 A, and L = 800 ^m-^ 2000 /im. According to Ref. [33], the frequency exponent for 
samples on a silicon substrate exceeds noticeably that for equally sized samples on a sapphire substrate (the ratio 
of (5's for the two cases is about 1.5). It follows then from the above formulas that the slope of the noise amplitude 
considered as a function of the film thickness is larger for samples on the silicon substrate. This is indeed the case as 
is clearly seen from Fig. 3 of Ref. [13] ■ Furthermore, since 6 increases for decreasing a, the slope of the noise amplitude 
considered as a function of the sample length is expected to be larger for thicker samples, but with a less noticeable 
difference in the slopes since the length exponent (6—1) is less sensitive to variations in 6 than the thickness exponent 
(—25) ((5 is about some tenths in both cases). This is again in agreement with the observations as is evident from 
Fig. 4 of Ref. [32] where the amplitude curve for a 800 A-thick film is somewhat steeper than that for a 400 A-thick 
film (both on a silicon substrate). 

Next, it is generally agreed that, with other things being equal, flicker noise is more intensive in semiconductors than 
in metals, and this is again in conformity with Eq. (|25p . because charge carrier mobility is higher in semiconductors 
than in metals, usually by several orders. Unfortunately, determination of mobility in semiconductors (or semimetals) 
is a difficult problem, both theoretically and experimentally, and different experiments often give significantly different 
results. By this reason, the subsequent consideration will be carried out for metals only. Even in this case careful 
estimation of the noise level takes some effort. This is because electron mobilities in thin metal films commonly used 
in fiicker noise measurements differ essentially from the corresponding bulk values, varying non-monotonically with 
the film thickness, and exhibiting complicated temperature dependence. Thus, the thicker the film, the more reliable 
comparison of theoretical and experimental results. Fortunately, the modern instrumentation allows measurements in 
sufficiently thick samples, electrical transport in which has bulk properties (usually, effects related to film thickness 
become important for h less than a few hundred nanometers). As is well known, temperature dependence of the 
electron mobility in this case is well approximated by the 1/T law. Theoretically, this approximation is valid for T 
higher than the Debay characteristic temperature, but in most cases it is practically applicable already for T > 50°K. 
Furthermore, the effective dimensionality D 3 in thick samples. Therefore, it follows from Eq. (p5| that in 
sufficiently thick samples the frequency exponent 7 = 1 + 5 ~ 1, and the flicker noise level is temperature independent. 
This conclusion is confirmed, e.g., by the results of Ref. [s^] where 1// noise was measured in 2.44 //m thick metal 
films, which is quite sufficient for bulk treatment of the sample conduction. According to Fig. 5 of Ref. 341, the fficker 
noise level is constant for T > 50°K indeed, and 7 is found to be about 1.01 for most samples.'^ In the opposite case 
of extremely thin films, conductivity is approximately independent of temperature. For instance, resistivity of a 5 nm 
thick gold film varies from the value 2.85 • lO'^Q • m at T = 50 °K to 3.0 • lO'^n • m at T = 275 °K, i.e., only by about 
5% ^36j . If (5 = 3 — 2) were constant, this would mean that the noise magnitude is a linear function of temperature in 



^ As mentioned in Ref. [sSl . at very high current densities excee ding 1.1 • lO'' A/cm^, the samples undergone structural defects resulting 
in somewhat higher values of 7. Unfortunately, the authors of !34fl did not specify the metals used in their experiments, which makes 
quantitative comparison with Eq. II25I I impossible. 
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this case. However, i5 is itself temperature dependent in thin films, namely, it is expected to decrease with increasing 
temperature, because thermal effects destroy the exchange correlations, thus raising the effective dimensionality of 
the sample. 

Finally, it was found in ^32"| that the frequency exponent in some cases depends also on the sample length, though 
much more weakly than on the sample thickness. This dependence cannot be explained from the point of view of the 
developed theory, so its appearance can serve as an indication on the limits of applicability of the theory. 



2. Quantitative comparison with the experiment 

In order to compare the absolute value of the noise spectral density given by Eq. ([25ljwith experimental data, we 
refer to the results of flicker noise measurements performed by Wong, Cheng and Ruan 33] , which were used already 
in the above qualitative analysis, and by Voss and Clarke Before going into detailed comparison, let us mention 



the following important circumstance. As is seen from Eq. ((25| . the noise amplitude is very sensitive to the value of 
S. This parameter appears, in particular, in the exponent of the ratio {c/ fa^) which is normally very large. Indeed, 
for / = 1 Hz and a = 10~^ cm, this ratio is equal to 3 • 10^°. Therefore, an error as small as 0.01 in the value of 5 
results in the extra factor of 10°'^ ~ 1.6 in the amplitude. At the same time, 6 is usually measured with the accuracy 
of a few hundredths at best, so the calculation carried out below is actually an order-of-magnitude estimation of the 
noise level. 

a. Wong, Cheng and Ruan. In this work, flicker noise power spectra were measured in the case of copper films 
of various thickness and length, sputtered on sapphire and oxidized silicon substrates. To illustrate the scheme of 
the calculation, let us take as an example the case of film with L — 1200 fim, a = 1,2- 10~^ cm, deposited on the 
sapphire wafer. According to Fig. 1 of [34|, samples of this thickness have conductivity a — 5.9 • 10® Q~^m~^ . Using 
the relation /i = a /en where n = 8 • 10^^ cm~'^ is the charge carrier concentration in copper, one finds the electron 
mobility fi = 1.4 • 10'^ units CGS. Next, according to Fig. 5 of [11], the frequency exponent 7 in the case under 
consideration is equal to 1.1, hence, (5 = 7—1 = 0.1. Substituting these values into the formula ([^5]) gives g = 510 
units CGS. Finally, putting this and** T = 300 °K in Eq. ([SO)) , one obtains k = 3.8-10"^^. This value is to be compared 
with the measured value >fexp — 6 • 10^^^ given in Fig. 3 of Ref. [12] (where it is called normalized noise amplitude). 
The experimental error of S is about 0.02, which implies an ambiguity by the factor of 2.5 in the calculated value of 

It was found in [12] that sufficiently thin copper films are characterized by a pronounced dependence of the frequency 
exponent on the sample length. Figure 5 of [33] shows that the thinner the film, the stronger this dependence: for 
L increasing from 800 /xm to 2000 /xm, 7 decreases by only 0.05 in the case of a = 800 A, and by 0.15 in the case of 
a — 400 A. As was already mentioned above, the present theory is unable to explain this dependence, and hence the 
corresponding experimental data are beyond the scope of applicability of the theory. By this reason, the quantitative 
comparison below is carried out only for films with a ^ 800 A. The calculated and measured values of k for films of 
various thickness on sapphire and silicon substrates are collected in Table I and Table II, respectively, together with 
the other parameters involved in the calculation. The quantities /x, (7, are given in the CGS system of units. 

It is seen from these tables that the calculated and measured values of k for films on sapphire wafer agree within 
the experimental error for all a. For films on silicon wafer, the agreement within the experimental error is found for 
a — 1200 Aand a = 1600 A. In the case of a = 800 A, where dependence of 7 on L is still noticeable, the theory 
somewhat overestimates the noise level. 

h. Voss and Clarke. In the work [3], fiicker noise was measured in thin metal films evaporated or sputtered on 
glass substrates. The information provided in this paper is sufficient for estimation of the noise intensity in the 
gold film shown in Fig. 2 of [2]. This was an elongated sample with h = 25 nm, w = 8/im, I = 625 /im, biased 
at Uq = 0.81V, and operated at about 40°K above room temperature. Unfortunately, the exact value of 7 is not 
given for this case by the authors, who mentioned only that it is close to 1. In view of what have been said about 
dependence of the noise amplitude on S, this implies a large amount of uncertainty in the theoretical estimation.^ 
Yet, assuming that ^ « indeed, we substitute the sample dimensions in Eq. and find g = 46cm~^. Next, 
in order to determine conductivity, we use the I — V characteristic of the given gold sample, shown in Fig. 3 of 
[2]. According to this figure, the sample resistance was about 100 f2. Taking into account the sample dimensions 



This value of T as well as the accuracy of 5, mentioned below, are communicated to the author by Prof. H. Wong. 
^ The measured spectrum presented in Q does not actually fit the l//'''-law even at low frequencies, because of the low-frequency roll 
off of the amplifier and capacitor used to improve impedance match for low-resistance samples. As a result, after subtraction of the 
background noise contribution, the low-frequency part of the corrected spectrum curve in Fig. 2 of Ref. Q goes above the measured 
values. 
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given above, this implies that a = 1.2 • 10^ ri~^ni~-'^. It should be mentioned that this value is approximately six 
times lower than that obtained in more recent studies of electrical transport in thin films. For instance, according to 
Ref. [s^l conductivity of a 25 nm thick, 15 fim wide gold film obtained by a laser-improved deposition of nanoparticle 
suspension, is 7.1 • 10^ n~^m~^. The same value can be obtained also indirectly using the data given in Refs. [3^,[3^. 
According to [3^, the conductivity of gold is 75% to 85% of its bulk value for h — 100 nm, depending on the choice of 
the substrate, and decreases below that value approximately linearly with decreasing thickness. On the other hand, 
according to Ref. [s^ conductivity drops to about 3 • 10^ il,~^m~^ for h = 5nm. One readily finds from this that for 
h = 25 nm, a = (6.5-^ 7.5) • 10^ il~^m~^. Presumably, this difference in the values of conductivity is to be attributed to 
the quality of film deposition. Substitution of cr = 1.2 • 10^ i7~^m~^ and n = 5.9 • 10^^ cm""^ in the relation fi — a /en 
yields the electron mobility ^ = 390 units COS. Then Eq. ([201) gives k ^ 9.6 ■ lO'^^ (for T = 330 °K). Putting this 
together with the bias value given above in Eq. (P5)l . we find \Cij\ — 6.3 • 10^^^ V^/Hz for the frequency / — 1 Hz, 
which is to be compared with the experimental value Cjj sa 10^^^ V^/Hz . 

V. DISCUSSION AND CONCLUSIONS 

We have shown that, staying within the one-particle picture of flicker noise generation by quantum electromagnetic 
fluctuations, backreaction of the conducting medium on the fluctuating field of the charge carrier can be described phe- 
nomenologically as an effective reduction of the system dimensionality. Using the dimensional continuation technique, 
we have found that the backreaction affects both the frequency dependence and the magnitude of the noise spectrum. 
Namely, the frequency exponent in the l//^-asymptotic of the fluctuation power spectrum in a conducting sample is 
found to be l-f S, where 6 — S — Ti being the effective dimensionality of momentum space, while dependence of the 
noise amplitude on S is given by Eqs. (j25p . (j26p . Although introduced initially as a momentum space characteristic, 
S thus relates the noise amplitude to geometric properties of the sample in the ordinary 3-dimensional coordinate 
space [Cf. Eqs. (p6|) . ((28)) ]. It was demonstrated in Sec. lIVBll that the experimentally observed dependence of the 
noise amplitude on the sample geometry for a given S is adequately described by Eq. ([^5]) . The way S itself depends 
on the sample thickness and system temperature, as expected from its definition, is also confirmed by observations. 
However, the value of S cannot be predicted within the developed approach. In other words, S plays the role of a 
phenomenological parameter of the theory. The noise amplitude turns out to be notably sensitive to the value of 
6: We saw in Sec. IIVB 21 that for a 10^^ cm-thick film, an increase of 0.1 in S raises the noise level by about two 
orders. This perfectly agrees with the observed rise of the noise level in samples with 7 > 1 in comparison^ with the 
predictions of the Hooge's empirical formula [s] obtained for 7 = 1. In fact, the quantitative comparison carried out 
in Sec. IIVB 21 shows that the calculated and measured values of the parameter x coincide within the experimental 
error. Deviations between the theory and experiment become noticeable only in films with a pronounced dependence 
of the frequency exponent on the sample length. It was mentioned in Sec. IIVB II that this dependence cannot be 
explained from the point of view of the developed theory, so its appearance indicates the limits of applicability of the 
theory. Although this issue is beyond the scope of the present approach, it is yet worth to comment on the possible 
origin of these deviations. According to Ref. [33 |. dependence of the frequency exponent on the sample length is 
noticeable in sufficiently thin films. At the same time, the fact that films deposited on different substrates exhibit 
different noise characteristics clearly shows that the boundary conditions affect significantly the mechanism of noise 
generation. From the theoretical point of view, this is reflected in the essential role played by the assumption that the 
film is plane-parallel in the derivation of the photon propagator [Cf. discussion below Eq. (jBSp j. On the other hand, 
this assumption is violated to some extent by the surface roughness of the film. In fact, the authors of 32J emphasize 
that the electric properties of the films used in their experiments are affected by the surface roughness, especially 
in the case of thin films. Thus, the surface roughness is a possible reason for the above-mentioned deviations in the 
values of >f. 

Finally, as was shown in Sec. IIV Al the obtained results explain the observed unboundedness of the flicker noise 
spectra, resolving naturally the problem of divergence of the total noise power. Together with the demonstrated qual- 
itative and quantitative agreement of the results with experimental data, this suggests that quantum electromagnetic 
fluctuations is the source of flicker noise in metal fllms. 



It is meant that compared are the noise levels at a fixed frequency, rather than the spectra themselves (in the latter case, of course, 
comparison would be meaningless). 
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APPENDIX A: GAUGE INDEPENDENCE OF THE POWER SPECTRUM 

Consider the theory of interacting scalar and electromagnetic fields described by the action 

S ^ + Sa , 

where is given by Eq. ([9]), and 

For arbitrary constant parameter a, the gauge fixing term describes the generalized Lorentz gauge. Let us introduce 
the generating functional of Green functions 

Z[J, T], jf] = J dAd(j>d(j>* exp + J d'^x[J''Af, + r;*0 + rjcj)*]^ | , (Al) 

where J,r],ri* denote sources for the fields A,(l)*,(p, respectively. Vanishing of Z under the gauge variation of the 
functional integral variables 

SAf, = df,^{x) , 6(j) = iei{x)(j) , 5dp* = -ie4{x)^* , (A2) 
with £^{x) a small gauge function, leads to the Ward identity 

Since we are interested in the connected contribution to the correlation function, we rewrite this identity for the 
generating functional of connected Green functions, W = —iln Z, 



The consequence of this equation we need is obtained by functional differentiation with respect to 77,77*, and twice 
with respect to J, with all the sources set equal to zero afterwards. 



SMy)SJa.{x)SMx')Si^iz)Sf]*{z') ' SJ^{x)Sj0{x')Sr]iz)Srj*{y) 



^eS^\y ^'^SMx)SJf3{x')dj^{y)dj^*{z') 
Fourier transform of this identity with respect to y reads 



SMy)SMx)SJ,3ix')Sr){z)Sr]*{z') 



/ —ik'z' —ik'z\ '"' ' • ( 

^ ~ 'sJa{x)SJp{x')Sr]{z)67j*{z') ' ^ ' 

The argument of the Fourier transform is purposely denoted here by k' to stress that the left hand side of this equation 
corresponds to the variation of the Green function we dealt with in Sec. IIIIl under gauge variation of the external 
field. Indeed, the longitudinal part of the photon propagator in the generalized Lorentz gauge has the form 

DUk) = -a^-^. (A6) 

Therefore, contraction with the factor k''^k'^/a is equivalent to amputation of the photon propagator attached to the 77 
vertex, followed by contraction of this vertex with k'^. Exactly the same result is obtained under the gauge variation of 
the external field coming into this vertex. The only difference with the Green function we considered in Sec. IIIIl is that 
the external scalar lines in Eq. (|A5p are the particle propagators. To promote them into particle amplitudes, according 
to the standard rules, Eq. (jASp is to be Fourier transformed with respect to the variables z, z' , and then multiplied by 
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ip{q)il^* {q'){m? — q^){m? — g'^), where the arguments q, q' of the Fourier transformations with respect to z, z' are to 
be taken eventuaUy on the mass sheU. But these operations give zero identicaUy when apphed to the right hand side 
of Eq. (jASp . because each of the factors e~*'^ ^ , e"*'^ ^ makes the corresponding particle propagator nonsingular on 
the mass shell. For instance, the first term in Eq. (jA5p gives rise to the contribution of the form {m? — q'^)D'^{q' + k') 
times terms nonsingular on the mass shell. For k' ^ 0, the function D'^{q' + k') is also nonsingular at g'^ = m^, and 
hence this contribution vanishes on the mass shell. 

Thus, the correlation function is invariant under gauge transformations of the external field, which are part of the 
gauge freedom in the theory. The other part is related to the explicit dependence of the photon propagator on the 
choice of the gauge conditions used to fix the gauge invariance of the action. As is well known, it is longitudinal 
part of the propagator that depends on the gauge. Let us first consider the simples case of Lorentz-invariant gauges. 
Then the most general form of the longitudinal part is given by Eq. (|A6[) in which a is to be regarded as an arbitrary 
function of k^ . It is not difficult to see that variations of a{k^) do not affect the observable quantities. Recall, first 
of all, that we are interested ultimately in the fluctuations of gauge-invariant quantities such as the electric field 
strength. The a-independence of these quantities is a direct consequence of their gauge invariance, because variations 
of a{k'^) give rise to terms that are pure gradients with respect to the spacetime arguments a;, a;', as is easily verified by 
substituting the expression (|A6p in place of one or two photon propagators in Eq. (|17p . Then, if the vector potential 
contribution to the field strength is negligible, as is the case in our nonrelativistic calculation (recall the condition 
|q| ^ m used throughout), the voltage correlation function can be found by integrating the correlation function for 
the field strength with respect to x^x' using the relation E — — VAq. 

More generally, the longitudinal part of the photon propagator in a Lorentz non-invariant gauge has the form, in 
coordinate space, 

where Xti,{^) is an arbitrary function of spacetime coordinates. If the Lorentz index of the spacetime derivative in 
dfj,Xv{^) is left free after combining the Feynman diagrams in Fig.[2l then this term leads to a gradient contribution to 
the two-point function of electromagnetic field, and, as before, does not contribute to the voltage correlation function. 
On the other hand, if the spacetime derivative is contracted with the interaction vertex, then the contribution of 
d^Xv{x) to the two-point function is not pure gradient, so the above argument does not work. Nevertheless, it can 
be shown that all such terms cancel each other in the complete expression for the correlation function. However, 
this requires examination of the complete set of Feynman diagrams, which complicates the proof. To avoid this 
complication, we will prove gauge independence of the low-frequency asymptotic of the power spectrum only, which 
is quite sufficient for our purposes. To this end, we use the following identity expressing invariance of the particle 
action under the transformation (IA2D 



5(t){x) ' S(t>*{x) " ' ' '"SA^ix) 
Differentiating this identity twice with respect to (piy), (t>*{z) and setting — afterwards yields 

^mTT^tIt^ttttt = ,./'ff , M ^\x -y)- f'^/,, M ^\x - z) . 
5A^{x)5(t){y)5(j)*{z) S(f>{x)6(f>* {z) ^ ^' Scj)* {x)S(f>{y) ^ ' 

The left hand side here is just the interaction vertex contracted with the derivative coming from the term 9^x^(a;). 
It follows that the result of this contraction is expressed through the inverse charged particle propagator. When built 
into the Feynman diagrams in Fig. [21 the latter is integrated with either the function il){z) (the external line), or the 
particle propagator. In the first case it gives zero by virtue of the definition of free particle states, while in the second 
it cancels the particle propagator. At the same time, as we saw in Sec. IIIIl it is singularity of the particle propagator 
that is responsible for the occurrence of l/w contribution to the function G(x^ a;', i', lS). Hence, longitudinal part of the 
photon propagator does not contribute to the leading low-frequency term of the power spectrum of electromagnetic 
fluctuations. 

Thus, gauge- independence of our results expressed by Eqs. ((24l) . ((25l) is proved completely. 



APPENDIX B: PHOTON PROPAGATOR IN THIN METAL FILM AT T = 0. 



For the sake of clarity, we keep track of the factor a in the formulas of the present section. As was mentioned 
in Sec. Ill Al the only component of the photon propagator, relevant in the calculation of the two-point correlation 
function, is Dqq{x — x') = i{'TAQ(x)Ao{x')). This important fact will be proved below in the case of the axial gauge 
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([2]) which was seen to be particularly convenient in describing the damping of the temporal photon propagation in 
the 3-dircction. Let the unit vector in this direction be denoted by n, 

To determine the form of the photon propagator at zero temperature 

D^,u{x - x')|t=o = D%{x - x') = {Q\7A^{x)A,{x'm , (Bl) 
we go over to momentum space 

D%(.^)^ I -0^e-^'^D%{k), (B2) 

and note that D'^j^{k) is a symmetric rank- two tensor that can be built only from the vectors n^, fc^, and the Minkowski 
tensor rj^^. Therefore, it has the following general structure 

where a, 6, c, and A are some scalars built from fc^, i.e., functions of the two invariant combinations [nk) = —k^ and 
k^ . According to discussion in Sec. Ill A[ the electromagnetic field operator is constrained by A3 = 0, dA^/dx^ = 0, 
which imply the following conditions on the form of the photon propagator 

r, I \ n dDQyjx) 

£'3^(a;) = 0, — = 0. (B3) 

0x3 

Written out in components, the first of these conditions reads 

{bko + ck3ko)A = 0, 
{bki + ck3ki)A = 0, « = 1,2, 
(-1 + a + 26fc3 + c/c|)A 0. 

It follows that 

a — 1 a — 1 



ck,^c{nk), c- 



^3 
so 



Then the second of the conditions (|B3[) gives 



{nk)r]o„ + (a - l)A;o"i. + (a - l)"^^^ | A = . 

The 3-component of this condition is satisfied identically, while the other three yield 

("fc) + («-l)(^}A = 0, (B5) 

(a-l)^A = 0, ^ = l,2. (B6) 

Since the photon propagator does depend on xo,xi,X2, it follows form Eq. (|B6[) that a = 1, and then Eq. (|B5[) gives 
k^A = implying that the function Z3^^(x) is independent of X3. Thus, 

D%{k)^iVf..+n^,n,)A. (B7) 
It remains only to find the quantity A. To this end, we will calculate the 00-component of the propagator explicitly. 
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Under the condition dAo/dx^ = 0, the normal mode decomposition of the scalar potential operator reads 
Ao{x, t) = N i — (flfe exp{i{-ujkt + km + k2X2)} + H.c.) , ujk = \k\ , 

where is a normalization factor to be determined below, and is the annihilation operator of a "temporal" photon 
with the wave vector k = (fci, fc2,0). As is known from quantum electrodynamics, the commutator of this operator 
with its Hermitian conjugate is afeOj^ — a^Cfc = — 1 . Substituting the expression for Ao in Eq. (|Bip . one finds 

DUx, t) = -iv2 ^ ^ exp{*(-a;fe|i| + kx)} . 

As always, the factor N depends on the choice of a "box" in which the field is quantized. Suppose that this box is a 
rectangular in (xi, a;2)-plane, its sides being parallel to the coordinate axes. For a sufficiently large box, summation 
over k can be replaced by integration over ScPk/{2Tr)'^, where S is the box area, 

-\-oo 

(Pk e*(-"'=l*l+'=^) f dko f (Pk e'^~''°*+'"^'> 



D°,{x,t) = -tN^S / — = N^S 



(27r)2 2|fe| J 2-K J (27r)2 fc2 - fc2 + iO ■ 

— OO 

A'^ can be found by calculating the electric potential produced by a resting point charge, e. If this charge is at the 
origin of the coordinate system, then the potential at the point x — (xi, X2, 0) is given by the well-known formula 

+ CXJ 

Ao{x) = ^e J dtD°„{x,t). 

— oc 

Substitution of the explicit expression for DQQ{x,t) yields 

/j2l. „ikx 

On the other hand, since this potential is independent of xa, the same result will be obtained if the film is replaced 
by an infinite set of identical parallel films adjacent to each other. This is allowed by our assumption that the film is 
plane-parallel. In other words, Aq(x) can be represented as a superposition of the Coulomb potentials produced by an 
infinite sequence of vertically aligned point charges, spaced at the distance a, or equivalently, by a charge distribution 
with density X2, X3) — {e/ a)5'''^\xi,X2)- Fourier transform of the latter is p{k) = [e/ a)2'K5{k^), hence 

. / N /",-?/ P{x') f 47re*'=^ 47re f d^k 



\x-x'\ J (27r)3'^' ' fc2 a J (27r)2 fc2 

and comparison with the preceding formula gives 

V oa 

The 00-component of the photon propagator thus takes the form 

'1*1' . 8»' «(*3) 



(27r)4 ' a k^ +iO 

Putting this A in Eqs. (|B7|) . (|B2p . we arrive finally at the following expression for the zero-temperature photon 
propagator 



I r d^k ,1.^ An 



e 



a J (27r)3 A;2 _^ jg 
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a, A 5_ /i X 10"^ G X 10"^ xth X lO" Xexp x 10 

800 0.14 1.1 9.8 15 9 

1200 0.10 1.4 5.1 3.8 6 

1600 0.10 2.2 4.7 5.5 5 



TABLE I: Calculated (xth) and measured (xexp) values of x for copper films of various thickness on sapphire substrates, /it, g, x 
are given in the CGS system of units. 



a, A r) // X iO a X iO - >^ai x il)-^- ^exp x 10 

800 0.21 0.9 32 22 2 

1200 0.17 1.2 14 4.9 2 

1600 0.14 2.0 8 2.3 1 



TABLE II: Same for copper films on silicon substrates. 
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FIG. 1: Symbolic diagrammatic picture of the effect of particle collisions and external electric field on the particle wave function. 
Wavy lines denote photon propagators, broken line the external electric field, solid lines charged particles. 
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FIG. 2: Feynman diagrams describing the first order external field correction to the particle propagator, q and p are the particle 
4-momentum and 4-momentum transfer, respectively. 



